Research Article 











Circulating angiogenic modulatory factors 
predict survival and functional class in 
pulmonary arterial hypertension 



Rajeev Malhotra 1 *, Samuel Paskin-Flerlage 2 *, Roham T. Zamanian 3 *, Patrick Zimmerman 4 , 
Jonathan W. Schmidt 4 , Donna Y. Deng 1 , Mark Southwood 5 , Robert Spencer 2 , Carol S. Lai 2 , William Parker 6 , 
Richard N. Channick 7 , Nicholas W. MorreW, C. Gregory Elliott 4 *, and Paul B. Yu 2 * 

department of Medicine, Cardiovascular Research Center, Cardiology Division, Massachusetts General Hospital, Harvard Medical 
School, 'Pulmonary and Critical Care Division, Massachusetts General Hospital, Harvard Medical School, Boston, Massachusetts, 
department of Medicine, Cardiovascular Division, Brigham and Women's Hospital, Harvard Medical School, Boston, Massachusetts, 
3 Vera Moulton Wall Center for Pulmonary Vascular Disease, Division of Pulmonary and Critical Care Medicine, Stanford University 
Medical Center, Stanford, California, "Department of Internal Medicine, University of Utah, Intermountain Medical Center, Murray, 
Utah, USA, department of Medicine, University of Cambridge School of Clinical Medicine, Addenbrooke's Hospital, Cambridge, UK 

and department of Surgery, Duke University Medical Center, Durham, North Carolina, USA 

'Authors contributed equally to this manuscript 



ABSTRACT 

The diagnosis of pulmonary arterial hypertension (PAH) is frequently delayed. We hypothesized that circulating angiogenic 
modulatory protein levels might correspond with vascular remodeling activity and serve as sensitive biomarkers of PAH. Levels 
of soluble endoglin (sEng), soluble vascular endothelial growth factor receptor- 1 (sVEGFRI), N-terminal brain natriuretic peptide 
(NT-proBNP), C-reactive protein (CRP), and other biomarkers were measured in peripheral blood from 97 PAH patients, 16 first- 
degree relatives of idiopathic or heritable pulmonary arterial hypertension (HPAH) patients, and 56 controls, and correlated with 
disease, functional class, hemodynamic parameters, exercise capacity, and transplant-free survival. Endoglin expression was 
analyzed in lung tissues of six individuals with idiopathic or HPAH and four individuals without PAH. Levels of sEng, sVEGFRI, 
CRP, and NT-proBNP were elevated in Group I PAH of diverse etiologies, with sEng performing better than NT-proBNP in detecting 
PAH (receiver operator characteristic-area-under-the curve [ROC-AUC] of 0.82 ± 0.03 vs. 0.71 ± 0.05, P = 0.016). While sEng, 
sVEGFRI, and NT-proBNP correlated with New York Heart Association (NYHA) class, sEng levels were more sensitive than 
NT-proBNP in detecting NYHA Class l-ll disease (ROC-AUC of 0.88 ± 0.05 vs. 0.67 ± 0.08, P = 0.028). sEng, sVEGFRI , CRP, 
and NT-proBNP predicted transplant-free survival by univariate Cox regression. After adjusting for NT-proBNP levels, each of the 
other three markers predicted transplant-free survival. In multivariate analysis, sEng and CRP were independent predictors of 
survival. Endoglin expression was markedly enhanced in the microvascular endothelium and endovascular lesions of PAH versus 
control lung tissues. Circulating angiogenic proteins sEng and sVEGFRI are sensitive markers of prognosis and function in Group 
I PAH, including mildly symptomatic disease, and may provide unique noninvasive data reflecting underlying remodeling activity. 
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Pulmonary arterial hypertension (PAH) is defined by an 
elevated mean pulmonary artery pressure (mPAP) of >25 
mmHg at rest with a left ventricular end diastolic pressure 
<15 mmHg, and a pulmonary vascular resistance (PVR] 
of >3 Wood units (WU). [1] PAH affects approximately 
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20 individuals per million adults with a 9-15% mortality 
rate in the first year despite current therapy 12 61 In the World 
Health Organization (WHO) classification, Group I PAH 
includes idiopathic pulmonary arterial hypertension (IPAH) 
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occurring in the absence of cardiopulmonary disease, 
a subset of which is heritable pulmonary arterial 
hypertension (HPAH], or associated pulmonary arterial 
hypertension (APAH) with connective tissue disease (CTD), 
stimulant drug use, congenital heart disease (CHD], portal 
hypertension, or human immunodeficiency virus (HIV]. [171 

Dysregulated endothelial cell and myofibroblast 
proliferation, smooth muscle hypertrophy, intimal fibrosis, 
and enhanced thrombosis — features evident in obstructive 
and obliterative arteriolar and multi-channeled "plexiform" 
lesions of PAH — have been interpreted by some to reflect 
an underlying process of disordered angiogenesis. [89] 
Endovascular PAH lesions are often accompanied by 
perivascular myelomonocytic infiltrates, suggesting a role 
of inflammation. 110121 While these features are found in 
end-stage PAH, it is unclear if dysregulated angiogenesis 
or vascular inflammation are present in early PAH or 
contribute to its progression. 

Recent observations have provided further evidence for 
dysregulated angiogenic activity in PAH. Administration 
of imatinib, an inhibitor of platelet-derived growth factor 
receptor beta (PDGF-R(B], which regulates smooth muscle 
proliferation and pericyte recruitment, attenuates PH in 
rats, and possibly in humans. [1315] In contrast, administering 
SU-5416, an antiangiogenic inhibitor of vascular endothelial 
growth factor receptor (VEGF-R] and PDGF-R(B signaling, 
triggers severe pulmonary vascular remodeling and 
pulmonary hypertension [PH] in rodents. 1161 The potent 
PDGF-Rp and SRC inhibitor dasatinib attenuates PH in 
rats, [17] whereas chronic dasatinib administration has 
been associated with the development of PAH in man. [18 ' 191 
Angiogenic signaling cascades have important, albeit 
complex, roles in the development and progression of PAH. 
We hypothesize that surrogates of angiogenic activity could 
provide insight into disease activity and progression as 
biomarkers of PAH. 

As noninvasive adjuncts to diagnostic right heart 
catheterization, circulating biomarkers have been evaluated 
in relationship to clinical severity, response to therapy, 
and survival in PAH, [1 ' 201 ideally providing additional risk 
stratification for clinical interventions. The most extensively 
validated PAH biomarkers, B-type natriuretic peptide 
(BNP] and its pro-peptide N-terminal brain natriuretic 
peptide (NT-proBNP], predict survival and function, but, 
as surrogates of right ventricular myocyte stretch, are 
insensitive for mildly symptomatic disease. 121231 Numerous 
other candidate biomarkers have been explored in PAH; 
however, their utility above and beyond BNP in prognosis 
and management have not been confirmed. [24] In a cohort 
of 97 individuals with Group I PAH and 56 age- and 
gender-matched unaffected controls, we evaluated several 
circulating markers of systemic and vascular inflammation, 



C-reactive protein (CRP], osteoprotegerin (OPG], and 
interleukin (IL]-15, and two markers of endothelial 
activation and angiogenic activity, soluble endoglin (sEng) 
and soluble VEGF receptor-1 (sVEGFRl]. We evaluated 
these biomarkers in relationship to functional status and 
survival, hypothesizing that surrogates of angiogenesis and 
inflammation might yield novel functional or prognostic 
data, and detect minimally symptomatic disease not 
captured by current biomarkers of PAH. 

MATERIALS AND METHODS 

Study population 

Plasma samples were obtained from 97 patients 
with confirmed diagnoses of PAH by WHO criteria 
(mPAP >2 5 mmHg at rest with a pulmonary capillary wedge 
pressure <15 mmHg and a PVR >3 WUs]. [11 This cohort 
included patients with idiopathic or HPAH, as well as PAH 
associated with CTD, CHD, portal hypertension, drugs 
and toxins, HIV, and Gaucher's disease. Samples were also 
collected from 16 first-degree relatives of IPAH or HPAH 
patients. Control samples were obtained from 56 healthy 
individuals. Patients, relatives, and healthy controls were 
recruited at the Pulmonary Hypertension Association 
Scientific Sessions, the California PH Forum Conference, 
and Stanford University Medical Center, with informed 
consent and Institutional Review Board (IRB] approval from 
Partners Health-care, Inc. (Boston, Mass., USA], University 
of Utah, LDS Hospital and Intermountain Medical Center 
(Salt Lake City, Utah, USA], and Stanford University Medical 
Center (Stanford, Calif., USA]. 

Serum samples and biomarker measurements 

Serum samples were collected from peripheral blood 
and stored at -80°C until analyzed. Serum analytes were 
measured by sandwich enzyme-linked immunosorbent 
assay (ELISA] according to manufacturers' protocols, 
including OPG (Duoset, R&D Systems, Minneapolis, Minn., 
USA], high-sensitivity CRP, sEng, sVEGFRl (Quantikine, 
R&D Systems], and IL-15 (Quantiglo, R&D Systems], or 
by competition ELISA in the case of NT-proBNP (Alpco 
Diagnostics, Salem, N.H.]. All analytes were measured at 
least in duplicate, and quantified in comparison to dilutions 
of recombinant protein standards and standardized serum 
samples. 

Clinical and hemodynamic parameters 

New York Heart Association (NYHA] functional class, 
6-Min Walk Distance (6MWD], pulmonary function 
tests, echocardiographic measurements of left and right 
ventricular function and Doppler-based estimates of 
right ventricular systolic pressure, medical therapies, and 
invasive hemodynamic data assessed contemporaneously 
with blood sampling (within three months before or after 
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sample collection] were obtained by chart review with 
informed consent. PVR was calculated as (mean pulmonary 
arterial pressure - pulmonary capillary wedge pressure]/ 
(cardiac output] and expressed in WUs. Survival data was 
obtained retrospectively both from hospital records and 
the social security death index. 

Immunohistochemistry 

Paraffin sections were obtained from explanted lung 
tissue of IPAH patients (n = 3), HPAH patients [n = 3], or 
controls (explanted lungs lacking pulmonary vascular 
disease, n = 2; untransplanted normal donor lungs, n = 1; 
and normal lung tissues, n = 1). After deparaffinization, 
antigen unmasking, and quenching of endogenous 
peroxidase activity (3% H 2 0 2 in methanol], slides were 
blocked with 5% normal goat serum in Tris-buffered 
saline with 0.1% Tween-20, incubated with biotinylated 
goat anti-human endoglin (R&D Systems, 1:50] overnight 
at 4°C, and incubated with streptavidin-horse radish 
peroxidase (LSAB2, Dako, Carpinteria, Calif, USA] at room 
temperature for 30 minutes followed by diaminobenzidine 
reagent (immPACT DAB, Vector Laboratories]. Sections 
were counterstained in hematoxylin, mounted, and imaged 
using an Olympus BX63 microscope equipped with a DP- 
25 color camera and cellSens Dimension 1.6 software 
(Olympus, Tokyo, Japan]. 

Statistical analyses 

Statistical analyses were performed using STATA 8.0 
(StataCorp, College Station, Tex., USA] or GraphPad 
Prism 5.0 (San Diego, Calif., USA]. Pairwise comparisons 
of plasma biomarker levels were made in the control 
population, PAH patients, and first-degree relatives using 
the Mann-Whitney statistic. Normality of biomarker levels 
and clinical variables were assessed by the Shapiro- 
Wilk test. Continuous variables were analyzed with 
either the Student's t-test or Mann-Whitney statistic 
as appropriate, and Fisher's exact test for categorical 
variables. Receiver operator characteristic (ROC] curves 
were derived and area-under-the curve (AUC] analysis 
performed to compare performance of biomarkers in 
predicting presence or absence of disease. Associations 
between biomarker levels were examined using 
Spearman's correlation. One-way analysis of variance 
(ANOVA] was used to determine associations between 
normally distributed (NT-proBNP, sEng, sVEGFRl, and 
OPG] biomarker levels and NYHA functional class, 
while Kruskal-Wallis testing was used for non-normally 
distributed biomarker levels (CRP and IL-15]. Survival 
characteristics were compared using the Kaplan-Meier 
method with Log Rank testing. Cox proportional hazards 
models were employed to identify univariate and 
multivariate predictors of transplant-free survival. Values 
are reported as mean + SEM unless noted, with P values 
<0.05 considered significant. 



RESULTS 

Study population 

Table 1 illustrates the baseline demographics of the 
control, PAH, and first-degree relative groups. The gender 
of patients in the PAH group reflected the known female 
predominance of PAH. Control individuals were selected to 
reflect a similar female predominance. The ages of the PAH 
and control populations were 48.9 ± 13.1 and 47.4 ± 8.0, 
respectively (P= 0.47). Among the disease group, 51.6% of 
cases were idiopathic (IPAH, n = 50], 21.6% associated with 
CTD (APAH-CTD, n = 21), 12.4% associated with stimulant 
drugs and toxins (APAH-Stim, n = 12], 6.2% associated 
with CHD (APAH-CHD, n = 6], and 8.2% associated with 
portal hypertension or other conditions (APAH-Other, 
n = 8). While the patients presented across the spectrum of 
functional status, the majority exhibited NYHA Class II or III 
symptoms at the time of evaluation (Table 2]. Hemodynamic 
parameters were obtained from right heart catheterization, 
with an average mPAP of 53 ± 13 mmHg, cardiac index 
of 2.7 ± 1.0 L/min/m 2 , and PVR of 11.5 ± 6.0 WUs in 
the 97 PAH patients. These patients received currently 
approved PAH-directed therapies including prostacyclin 
analogues, endothelin receptor antagonists (ETRA], and 
phosphodiesterase-type 5 (PDE5] inhibitors. More than 
50% received combination therapy and a similar proportion 
received chronic anticoagulation. 

Levels of CRP, NT-proBNP, sEng, and sVEGFRl, 
but not OPG or IL-15 were elevated in PAH 

Baseline levels of biomarkers in our control population fell 
within the range of normal values previously reported. 12125 281 
Serum levels of CRP (P< 0.0001), NT-proBNP (P= 0.0008), 
sEng (P< 0.0001), and sVEGFRl (P< 0.0001) were elevated 
in individuals with PAH compared to control individuals of 
similar age and gender (Fig. 1 and Table 3). However, no 
differences in the levels of OPG and IL-15 were observed 



Table 1 : Population demographics and pulmonary arterial 
hypertension etiology 





PAH 

patients 


Control 
patients 


1° Relatives of 
IPAH/HPAH 
Patients 


Total sample size, n 


97 


56 


16 


Gender, % female 


89.7 


86.3 


75.0 


Age, years±SD 


48.9±13.1 


47.4±8.0 


54.4±11.9 


PAH subtype, % (n) 








IPAH/HPAH 


51.6 (50) 






APAH-CTD 


21.6 (21) 






APAH-Stim 


12.4 (12) 






APAH-Other* 


8.2 (8) 






APAH-CHD 


6.2 (6) 







IPAH: idiopathic pulmonary arterial hypertension; HPAH: heritable 
pulmonary arterial hypertension; APAH: associated pulmonary arterial 
hypertension; CTD: connective tissue disease; Stim: stimulant drug and 
toxin; CHD: congenital heart disease; *Other includes PAH associated with 
portal hypertension (n=5); HIV (n=2), and Gaucher's disease (n = l); 
PAH: pulmonary arterial hypertension; HIV: human immunodeficiency virus 
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between the PAH cohort and controls. None of the PAH- 
directed therapies taken in isolation or in combination 
(i.e., prostacyclins, ETRA, or PDE5 inhibitors) were 
associated with biomarker levels, except that patients 



Table 2: Clinical, hemodynamic, and treatment 
characteristics of pulmonary arterial hypertension patients 
and the idiopathic pulmonary arterial hypertension 
subgroup 



Clinical variables 


All PAH 


IPAH/HPAH 




(17=97) 


(n=50) 


NYHA functional class* % (n) 






Class I 


5 (3) 


3 (1) 


Class II 


58 (35) 


60 (18) 


Class III 


28 (17) 


33 (10) 


Class IV 


8 (5) 


3 (1) 


6MWD (m) mean±SD 


424±101 


411+103 


Hemodynamic variables* mean±SD 






RAP (mmHg) 


11+5 


10+5 


PASP (mmHg) 


77±28 


78+35 


mPAP (mmHg) 


53±13 


52+13 


CO (L/min) 


4.3+1.5 


4.7+1.4 


Cardiac index (L/min/m J ) 


2.7+1.0 


3.0+1.0 


PVR (Wood units) 


11.5+6.0 


10.8+5.2 


PAH-directed treatment* % (n) 






Prostacyclin analogue 


50.7 (35) 


50.0 (16) 


Endothelin antagonist 


55.1 (38) 


50.0 (16) 


PDE-5 inhibitor 


49.3 (34) 


46.9 (15) 


Combination** 


56.5 (39) 


50.0 (16) 


None 


11.6 (8) 


15.6 (5) 



*NYHA functional class and treatment data were available for 60 and 69 of 
the 97 PAH patients, respectively, and for 30 and 32 of the 50 IPAH/HPAH 
patients; respectively. Six-min walk data within 3 months of blood sample 
collection were available in 51 of the PAH patients and 21 of the IPAH/HPAH 
patients. Hemodynamic values within 3 months of blood sample collection 
were available in 35 of the PAH patients and 18 of the IPAH/HPAH patients. 
**Receiving two or more PAH-directed therapies; NYHA: New York Heart 
Association; 6MWD: 6 min walk distance; PDE: phosphodiesterase; IPAH: 
idiopathic pulmonary arterial hypertension; HPAH: heritable pulmonary 
arterial hypertension; APAH: associated pulmonary arterial hypertension; 
PAH: pulmonary arterial hypertension; PVR: pulmonary vascular resistance; 
RAP: right atrial pressure; mPAP: mean pulmonary artery pressure; 
CO: cardiac output; PASP: pulmonary artery systolic pressure 



receiving prostacyclin analogues exhibited reduced 
sVEGFRl levels (121 ± 10 vs. 172 ± 14 pg/mL, P = 0.005). 

Levels of sEng were more sensitive and specific 
for PAH than NT-proBNP 

Consistent with prior reports, 121231 elevated NT-proBNP 
predicted PAH with an ROC-AUC of 0.71 ± 0.05 (Fig. 2). 
Only sEng was superior to NT-proBNP in predicting the 
presence of PAH (AUC 0.82 ± 0.03, P = 0.016), while 
sVEGFRl (0.79 ± 0.05) and CRP (0.77 ± 0.04) were similar 
to NT-proBNP. An equally weighted combination of sEng and 
NT-proBNP (calculated as the linear sum of each biomarker 
value's Z score, i.e., [value - mean] /(standard deviation) 
yielded an AUC of 0.87 + 0.04, which was not significantly 
better than sEng alone. 

Relationship of biomarker levels and underlying 
etiology of PAH 

This cohort included patients with distinct etiologies 
of WHO Group I PAH disease (Table 1). The mean and 
distribution of each biomarker studied in relationship to 
underlying disease etiology is shown numerically (Table 3) 
and graphically (Fig. 3). Circulating levels of sVEGFRl were 
elevated in all categories of Group I PAH, whereas levels 
of sEng and CRP were increased in all groups studied 
(IPAH, APAH-CTD, APAH-toxin, and APAH-other) with the 
exception of CHD. As observed among PAH patients as a 
whole, among IPAH and HPAH patients sEng performed 
better than NT-proBNP in predicting disease (AUC of 
0.86 ± 0.04 vs. 0.63 ± 0.07, P = 0.008). 

Levels of OPG and IL-15 were elevated in 
APAH-CTD 

Elevated OPG levels were not associated with PAH; 



Table 3: Serum levels of inflammatory markers, angiogenic modulators, and natriuretic peptides in controls, pulmonary 
arterial hypertension patients, and 1 ° relatives of idiopathic pulmonary arterial hypertension 



Group 


CRP 
(ng/mL) 


NT-proBNP 
(fmol/mL) 


sEng 
(ng/mL) 


sVEGFRl 
(pg/mL) 


OPG 
(pg/mL) 


IL-15 
(pg/mL) 


Control 


713+131 (54) 


155.7+7.3 (28) 


3.66+0.12 (56) 


98.7+5.6 (39) 


2891+156 (55) 


1.88+0.13 (28) 


1° Relatives 


810+182 


174.4+32.6 


3.75+0.26 


121.1+10.4 


2525+413 


2.48+0.74 


of IPAH/HPAH 


P=0.09 (16) 


P=0.78 (11) 


P=0.72 (16) 


P=0.15 (6) 


P=0.11 (16) 


P=0.99 (11) 


All PAH 


2208+308 


229.3+14.4 


5.15+0.15 


141.6+5.98 


3392+156 


2.15+0.21 




P<0.0001 (97) 


P=0.0008 (76) 


P<0.0001 (97) 


P<0.0001 (68) 


P=0.09 (95) 


P=0.64 (42) 


IPAH/HPAH 


2066+397 


196.0+14.8 


5.12+0.19 


132.1+6.3 


3211+189 


1.68+0.16 




P<0.0001 (50) 


P=0.08 (36) 


P<0.0001 (50) 


P=0.0003 (33) 


P=0.31 (49) 


P=0.17 (17) 


APAH-CTD 


2181+553 


264.5+38.7 


5.43+0.40 


135.8+14.7 


3917+416 


3.60+0.87 




P=0.0007 (21) 


P=0.0009 (14) 


P<0.0001 (21) 


P=0.0047 (17) 


P=0.0239 (21) 


P=0.0023 (8) 


APAH-Stim 


3112+1440 


330.5+53.9 


5.26+0.38 


194.0+20.8 


3232+506 


1.83+0.17 




P=0.0013 (12) 


P=0.0002 (12) 


P=0.0001 (12) 


P=0.0002 (6) 


P=0.72 (11) 


P=1.0 (10) 


APAH-CHD 


1368+729 


160.0+21.5 


3.95+0.51 


145.3+13.5 


3325+546 


1.23+0.41 




P=0.28 (6) 


P=0.75 (6) 


P=0.89 (6) 


P=0.0084 (4) 


P=0.44 (6) 


P=0.17 (3) 


APAH-Other 


2437+1092 


217.3+21.4 


5.30+0.33 


152.4+21.8 


3393+514 


2.70+0.26 




P=0.0046 (8) 


P=0.0087 (8) 


P=0.0001 (8) 


P=0.0039 (8) 


P=0.45 (8) 


P=0.0263 (4) 



Values are presented as mean+SEM (n values in parentheses). All P values were determined using the Mann-Whitney statistic comparing the disease or first-degree 
relative group with the control group for that particular biomarker; OPG: osteoprotegerin; NT-proBNP: N-terminal brain natriuretic peptide; sEng: soluble 
endoglin; sVEGFRl: soluble VEGF receptor-1; CRP: C-reactive protein; IPAH: idiopathic pulmonary arterial hypertension; HPAH: heritable pulmonary arterial 
hypertension; APAH: associated pulmonary arterial hypertension; PAH: pulmonary arterial hypertension; CTD: connective tissue disease; Stim: stimulant drug 
and toxin; CHD: congenital heart disease 
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Figure 1: C-reactive protein (CRP), N-terminal brain natriuretic peptide (NT-proBNP), soluble endoglin (sEng), and soluble VEGF receptor- 1 (sVEGFRl) are 
elevated in pulmonary arterial hypertension (PAH). Levels of the biomarker analyzed are depicted among control individuals, PAH patients, and the first-degree 
relatives of either idiopathic pulmonary arterial hypertension or heritable pulmonary arterial hypertension patients, with median, top and bottom quartile values 
indicated as shown. 
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Figure 2: Soluble endoglin (sEng) is a better predictor of the presence of pulmonary arterial hypertension than N-terminal brain natriuretic peptide (NT-proBNP). 
Receiver operator characteristic curves for sEng, soluble VEGF receptor- 1 (sVEGFRl), C-reactive protein (CRP), and osteoprotegerin (OPG) are depicted in 
relationship to NT-proBNP as a reference biomarker. Areas under the curve (AUC ± SEM) were determined for sEng (0.824 ± 0.033), sVEGFRl (0.785 ± 0.045), 
CRP (0.769 ± 0.042), NT-proBNP (0.714 ± 0.051), and OPG (0.584 ± 0.049), T<0.05 versus NT-proBNP; # P<0.05 versus OPG. 
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Figure 3 : Angiogenic biomarkers are elevated in diverse etiologies of group I pulmonary arterial hypertension (PAH). Levels of inflammatory markers (C-reactive 
protein [CRP], IL-15, and osteoprotegerin [OPG]), angiogenic markers (soluble endoglin [sEng] and soluble VEGF receptor- 1 [sVEGFRl]) and N-terminal 
brain natriuretic peptide (NT-proBNP) levels are depicted among the control population, and patients with PAH of various etiologies: idiopathic or heritable 
PAH (IPAH), patients with PAH associated with connective tissue disease (CTD), stimulant drug or toxin (Stim), congenital heart disease (CHD), or other forms 
(portopulmonary, human immunodeficiency virus, and Gaucher's disease) of PAH (*P<0.05, f P<0.01, *P<0.001 as compared to control group), with median, 
top and bottom quartile values indicated as shown. 



however, among APAH-CTD patients, OPG levels were 
increased (3917 ± 416 pg/mL vs. 2891 ± 156 pg/mL in 
controls, P = 0.02). Similarly, elevated IL-15 levels were 
not associated with PAH, but were increased in APAH-CTD 
patients (3.60 ± 0.87 vs. 1.88 ± 0.13 pg/mL, P = 0.002], 
underscoring the known association between IL-15 and 
rheumatologic disorders (Fig. 1 and Table 3). [27 ' 281 By 
ROC analysis, OPG and IL-15 levels were sensitive for the 
presence of PAH among individuals with APAH-CTD (AUCs 
of 0.67 ± 0.07 and 0.86 ± 0.07, respectively). 

Biomarkers were not elevated in first-degree 
relatives of patients with IPAH or HPAH 

The levels of each of the biomarkers tested were not 
different between the first-degree relatives of IPAH and 
HPAH patients and controls (Table 3). 

Association of biomarker levels with functional 
class 

Elevated NT-proBNP (P = 0.03), sEng (P = 0.02), and 
sVEGFRl (P= 0.004) levels were associated with increasing 
NYHA Class, while OPG (Fig. 4), CRP, and IL-15 levels were 
not (data not shown). Levels of NT-proBNP, sEng, and 
sVEGFRl were elevated among patients of NYHA Class I or II 
status; however, sEng was a significantly better predictor 



of the presence of disease among NYHA Class I - II patients 
than NT-proBNP by ROC analysis (AUC 0.88 ± 0.05 vs. 
0.67 ± 0.08 respectively, P = 0.028). Among these markers, 
NT-proBNP levels also demonstrated an association with 
6MWD (r= -0.37, P= 0.018). 

sEng correlates highly with sVEGFRl, but not 
NT-proBNP and CRP 

Among individuals with PAH, the strongest correlation 
between the biomarkers existed between levels of sEng and 
sVEGFRl (r = 0.39, P = 0.001, Table 4, Fig. 5). In contrast, 
sEng did not correlate with CRP or NT-proBNP. 

Predictors of survival in PAH 

Median follow-up for the PAH patients was 3.2 years 
(interquartile range, IQR 3.0-4.2), while 5-year transplant- 
free survival was 73%. Of clinical and hemodynamic 
variables, univariate predictors of transplant-free survival 
included NYHA, 6M WD, right atrial pressure, cardiac index, 
and PVR (Table 5). Age was not a predictor of transplant- 
free survival in our PAH population (Cox HR 1.03 per 
decade, P = 0.88). Analyzed as continuous variables, 
biomarkers that predicted transplant-free survival included 
NT-proBNP, sEng, sVEGFRl, and CRP with Cox hazard 
ratios (HR) of 1.43 (for every 100 fmol/mL increase, 95% 
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Figure 4: Soluble endoglin (sENG), soluble VEGF receptor- 1 (sVEGFRl), and N-terminal brain natriuretic peptide (NT-proBNP) levels in pulmonary arterial 
hypertension (PAH) patients are associated with New York Heart Association (NYHA) Class. Levels of biomarkers are depicted for control individuals and PAH 
patients grouped by NYHA Class, with median, top and bottom quartile values indicated as shown. Trend analysis was performed with the one-way analysis 
of variance for the diseased population. 
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Figure 5: Soluble endoglin (sEng) levels correlate with levels of soluble VEGF receptor- 1 (sVEGFRl ). (A-C) Circulating levels of sEng correlated with levels 
of sVEGFRl (r=0.39, i*=0.001), but did not correlate with levels of C-reactive protein (CRP) or N-terminal brain natriuretic peptide (NT-proBNP). 



confidence interval [CI] 1.11-1.85, P = 0.005], 1.42 (for 
every 1.0 ng/mL increase, 1.08-1.87, P = 0.013], 1.14 (for 
every 10 pg/mL increase, 1.02-1.27, P = 0.025], and 1.11 
(for every 500 ng/mL increase, 1.05-1.17, P < 0.001], 
respectively. Each of these biomarkers remained as 
significant predictors of transplant-free survival after 
adjusting for age and gender. 

Kaplan-Meier survival analysis was performed on 
biomarkers dichotomized about the median value 
(Fig. 6]. Median values for NT-proBNP, sEng, sVEGFRl, 



and CRP were 200 fmol/mL, 4.97 ng/mL, 131 pg/mL, 
and 1100 ng/mL, respectively. PAH patients with an 
NT-proBNP <200 fmol/mL exhibited a 74% reduced 
hazard rate for death or transplant (Cox HR 0.26, 95% CI 
0.08-0.82, P = 0.02] compared to those with NT-proBNP 
>200 fmol/mL. Similarly, PAH patients with either a 
sEng level <4.97 ng/mL, sVEGFRl <131 pg/mL, or CRP 
level <1100 ng/mL exhibited improved transplant-free 
survival with a Cox HR of 0.33 (95% CI 0.11-1.0, P = 0.05], 
0.13 (95% CI 0.02-1.0, P = 0.05], and 0.22 (95% CI 0.06- 
0.78, P= 0.02], respectively. 
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Table 4: Correlations among pulmonary arterial 
hypertension biomarkers 



Biomarker 



CRP 



NT-proBNP 



sEng 



sVEGFRl 



OPG 

sVEGFRl 
sEng 

NT-proBNP 



-0.10 

P=0.35 
(95) 
0.07 

P=0.58 
(68) 
0.14 

P=0.18 
(97) 

0.22 

P=0.06 

(76) 



0.14 

P=0.22 
(74) 
0.38 

P=0.005 
(55) 
0.19 

P=0.10 
(76) 



0.04 
P=0.68 

(95) 

0.39 
P=0.001 

(68) 



0.05 
P=0.70 
(67) 



Spearman's rho correlation with associated P values; sample size indicated 
in parentheses; OPG: osteoprotegerin; NT-proBNP: N-terminal brain 
natriuretic peptide; sEng: soluble endoglin; sVEGFRl: soluble VEGF 
receptor-1; CRP: C-reactive protein 



Table 5: Univariate predictors of transplant-free survival 
in pulmonary arterial hypertension 



Variable 


Cox hazard ratio 


P value 




(95% confidence 






interval) 




NYHA (for each class increase) 


4.71 (2.17-10.25) 


<0.001 


6 min walk distance 


1.56 (1.00-2.45) 


0.05 


(for every 50 m 1) 






Right atrial pressure 


2.57 (1.32-5.01) 


0.006 


(for every 5 mmHg T) 






Cardiac index 


4.45 (1.88-10.57) 


0.001 


(for every 0.5 L/min/m 2 i) 






Pulmonary vase, resistance 


1.16 (1.05-1.29) 


0.004 


(every 1 WU T) 






NT-proBNP 


1.43 (1.11-1.85) 


0.005 


(every 100 fmol/mL T) 






sEng (every 1.0 ng/mL T) 


1.42 (1.08-1.87) 


0.013 


sVEGFRl (every 10 pg/mL T) 


1.14 (1.02-1.27) 


0.025 


CRP (every 500 ng/mL T) 


1.11 (1.05-1.17) 


<0.001 



NYHA: New York Heart Association; WU: wood unit; NT-proBNP: N-terminal 
brain natriuretic peptide; sEng: soluble endoglin; sVEGFRl: soluble VEGF 
receptor-1; CRP: C-reactive protein 



In multivariate Cox regression analysis, sEng and CRP 
were independent predictors of transplant free survival 
with a Cox HR of 1.40 (for every 1.0 ng/mL increase, 1.04- 
1.90, P = 0.028) and 1.10 (for every 500 ng/mL increase, 
1.04-1.16, P < 0.001], respectively. Right atrial pressure 
(RAP) correlated with NT-proBNP (r = 0.42, P = 0.02) 
but did not correlate with any of the other biomarkers 
suggesting that sEng, sVEGFRl, and CRP predict prognosis 
in PAH independent of hemodynamic right ventricular 
compromise. Furthermore, after adjusting for NT-proBNP 
levels, sEng, sVEGFRl, and CRP each predicted transplant- 
free survival with Cox regression analysis. 

Endoglin expression is enhanced in human 
endovascular lesions of PAH 

Given that sEng may derive from the cleavage of membrane 
endoglin in endothelial cells, 1291 we examined expression of 
endoglin in lung tissues from patients with IPAH or HPAH, and 



II 



individuals without pulmonary vascular disease. Endoglin 
was expressed prominently in the endothelium of large, 
medium, and small arteries, arterioles, and venules, and to 
a much lesser extent in the adventitia, but not the media of 
the pulmonary vessels (Fig. 7A). No difference in endoglin 
staining between patients and controls was observed in 
large and medium-sized pulmonary arteries, whereas small 
pulmonary arteries (diameter < 40 \im) and peri-bronchiolar 
arteries exhibited increased endoglin expression in IPAH 
and HPAH (Fig. 7B and C). Most impressively, intraluminal 
channels of plexiform lesions exhibited intense expression 
of endoglin in both IPAH and HPAH lung tissues. 

DISCUSSION 

PAH is a multifactorial disease involving abnormal vascular 
tone, endothelial dysfunction, inflammation, dysregulated 
angiogenesis, and enhanced thrombosis. 19101 In this study, 
we postulated that circulating markers of angiogenic 
activity or inflammation might serve as sensitive markers of 
PAH disease severity and prognosis, and due to their ability 
to reflect vascular remodeling activity, perform better than 
conventional biomarkers. Widely used markers such as BNP 
and NT-proBNP correlate closely with functional status, 
therapeutic responses, and outcomes in PAH, [22;231 but are 
limited by the fact that right ventricular distention reflected 
by these markers occurs in relatively advanced disease. 
In fact, the potent anti-angiogenic molecule sEng had the 
best overall performance based upon ROC analysis among 
mildly symptomatic (NYHA Class I-II) patients, suggesting 
sensitivity even with modest disease burden. sEng, 
sVEGFRl, and NT-proBNP predicted functional class and 
disease severity and, along with CRP, predicted transplant- 
free survival in PAH by univariate analysis. After adjusting 
for NT-proBNP levels, sEng, sVEGFRl, and CRP each still 
predicted survival in PAH, indicating that these biomarkers 
predict prognosis in PAH above and beyond right 
ventricular compromise. Only sEng and CRP independently 
predicted survival based on multivariate analysis, however. 
Importantly, endoglin expression was enhanced in the 
microvascular endothelium of IPAH and HPAH lung tissues 
as compared to non-PAH control lungs, and was intensely 
expressed in plexiform lesions. Endothelial-derived anti- 
angiogenic proteins such as sEng and sVEGFRl, known 
to be released in the contexts of endothelial activation 
and vascular remodeling, appear to be more sensitive and 
specific in Group I PAH than previously described markers, 
and may represent surrogates of angiogenic remodeling 
activity in PAH. The elevation of anti-angiogenic markers 
in Group I PAH of diverse etiologies suggests dysregulated 
angiogenesis may be a common feature of these conditions. 

Endoglin and VEGFR1 are expressed in the vascular 
endothelium where they transduce prosurvival and 
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Figure 6: Soluble endoglin (sEng), soluble VEGF receptor- 1 (sVEGFRl), N-terminal brain natriuretic peptide (NT-proBNP), and C-reactive protein (CRP) 
predict transplant-free survival in pulmonary arterial hypertension. Kaplan-Meier survival curves are depicted for sEng, sVEGFRl, NT-proBNP, and CRP, all 
of which were dichotomized about their respective median value in our cohort. Median values for sEng, sVEGFRl, NT-proBNP, and CRP in this pulmonary 
arterial hypertension population were 4.97 ng/mL, 131 pg/mL, 200 fmol/mL, and 1100 ng/mL, respectively. Transplant-free survival curves were compared 
using the Log Rank test. 



proproliferative signals required for vascular homeostasis 
and angiogenesis. [30 32] In contrast, the soluble extracellular 
domains of endoglin and VEGFR1 appear to exert potent 
antiangiogenic effects. Endoglin, induced with endothelial 
activation, serves as a coreceptor for bone morphogenetic 
protein (BMP] ligands such as BMP9 and BMP10 which 
promote endothelial survival and quiescence. 133341 sEng, 
which blocks signaling of these ligands, potently inhibits 
tumor xenograft vascularization in vivo. [35] Similarly, 
sVEGFRl is a truncated splice variant of VEGFR1 that 
inhibits angiogenesis and induces endothelial dysfunction 
by sequestering VEGF and placental growth factor. [32,361 
Both sEng and sVEGFRl are markedly elevated in patients 
with pre-ecclampsia, and in addition to being sensitive 
biomarkers of this disease, appear to contribute to the 
underlying pathophysiology. 125 371 The elevation of sEng and 
sVEGFRl seen in PAH may also reflect endothelial activation 
and proliferation common to these two conditions. While 



sVEGFRl has been found to be elevated in patients with PH 
associated with sickle cell disease, 1381 and sEng was found to 
be elevated among individuals with scleroderma and PAH, we 
found increased levels of sEng and sVEGFRl in diverse Group 
I PAH patients, including those with minimal symptoms. 
The observation of markedly enhanced endoglin expression 
in plexiform lesions, thought to be sites of dysregulated 
angiogenesis in PAH, [9101 would support the interpretation 
that circulating sEng reflects endothelial activation, 
proliferation, and active pulmonary vascular remodeling. 
sEng may be released from sites of active remodeling by the 
activity of matrix metalloproteinase- 14 (MMP-14], reported 
to regulate endoglin cleavage in carcinoma tissues. 1291 We did 
not find elevated vascular expression of MMP-14 in normal 
or PAH lung tissues (data not shown], suggesting that other 
mechanisms might contribute. Endoglin is expressed in the 
failing left ventricle, and its soluble isoform is detected in 
the circulation after subjecting the left ventricle to pressure 
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Figure 7: Endoglin expression is increased 
in the pulmonary vascular endothelium of 
pulmonary arterial hypertension (PAH) 
patients. (A) Immunofluorescent staining 
of normal human lung reveals endoglin 
expression (green) in the endothelium of 
muscularized arterioles (ma) and non-MA 
and venules (nm), which does not overlap 
with expression of vascular or bronchial 
(br) smooth muscle a-actin depicted in red 
(scale bar=50 urn). (B) Endoglin (stained in 
brown using 3,3'-diaminobenzidine, DAB, 
with hematoxylin counterstain in blue) 
was expressed in the endothelium of large, 
medium and small vessels of control lungs, 
expressed to an equal extent in the large and 
medium vessels of idiopathic pulmonary 
arterial hypertension (IPAH) and heritable 
pulmonary arterial hypertension (HPAH) 
lungs, enhanced in the small (diameter 
<40 um) vessels of IPAH and HPAH lungs, 
and markedly expressed in the luminal 
endothelium of plexiform lesions (scale 
bar=40 um). (C) Endoglin expression, 
quantified by DAB staining intensity, was 
substantially increased in vessels with a 
diameter of <40 um in lungs of individuals 
with IPAH or HPAH, as compared to vessels 
of the same size analyzed in lung tissues 
obtained from individuals with airway 
disease without PAH (values shown as 
mean±SEM, PO.001 vs. controls). 



overload. 1391 Circulating endoglin can also be detected in 
the context of endothelial microparticles, as others have 
observed in PAH. [401 Our current approach, however, did not 
distinguish between the expression of sEng as a secreted 
isoform versus proteolytic cleavage product, or as a 
ventricular- or endothelium-derived fragment. 

In our study, levels of circulating sEng and sVEGFRl in 
PAH patients were 2- to 3-fold increased as compared to 
controls, lower than levels reported in pre-ecclampsia. [2SI 
If these antiangiogenic proteins derive from a pulmonary 
vascular bed that is obliterated due to extensive remodeling 
in PAH, the observed increases in these proteins may reflect 
proportionally greater expression by a diminished arteriolar 
tree. Since hypoxia itself can induce endoglin cleavage, 
elevated levels of endoglin associated with PAH could simply 
reflect hypoxemia. We did not observe increased expression 
of endoglin in the vessels of lungs obtained from individuals 
with known hypoxemic airway disease without PAH, as 
compared to lung tissues without airway disease or PAH, 
decreasing the chance of such a nonspecific finding (Fig. 7). 

The performance of sEng as a biomarker might suggest a 
pathophysiologic role in the progression of PAH, but further 
work in appropriate models will be needed to ascertain 



whether antiangiogenic factors such as sEng are causal, 
protective, or secondary to pulmonary vascular remodeling. 
Future studies are also needed to determine whether risk 
stratification based on the novel antiangiogenic biomarkers 
should impact the therapeutic approach taken for individual 
patients. 

Elevated OPG levels have previously been reported 
among individuals with IPAH, [41421 whereas OPG was 
only elevated among APAH-CTD patients in our cohort. 
OPG levels predict mortality associated with systemic 
cardiovascular disease, atherosclerosis, or diabetes 
mellitus. [26 ' 43 46] OPG, being a key modulator of tumor 
necrosis factor-a signaling and osteoclast function, 
may reflect states of vascular inflammation and bone 
remodeling, as well as vascular calcification, itself a 
marker of cardiovascular risk. [47 ' 48] IL-15, an inflammatory 
cytokine observed to be elevated in systemic sclerosis 
with thoracic involvement, 1271 is also thought to contribute 
to the pathophysiology of atherosclerosis and systemic 
vascular disease. [49 51] Since the levels of these markers 
were not measured in CTD patients lacking PAH, it is 
possible that the elevation of OPG and IL-15 among APAH- 
CTD patients reflects systemic inflammation associated 
with autoimmunity. 
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In contrast, CRP, an acute phase reactant and established 
biomarker of cardiovascular risk used in stratification of 
patients for coronary disease primary prevention, [52] was 
a sensitive marker of Group I PAH disease, consistent with 
previous reports. 1531 While some groups have not found CRP 
to predict survival in PAH, [541 our study validates CRP as an 
important predictor of long-term prognosis in PAH. CRP 
and sEng were independent predictors of transplant-free 
survival, suggesting that these markers may reflect distinct 
aspects of PAH progression or severity, a notion supported by 
the lack of correlation between these two markers (Table 4). 

None of the biomarkers assayed were significantly elevated 
in the relatively small cohort of asymptomatic first-degree 
relatives of IPAH and HPAH patients, suggesting these 
biomarkers may be surrogates of disease activity rather 
than disease potential. However, several of the biomarkers 
exhibited variability among controls and first-degree 
relatives, with CRP trending toward being elevated in first- 
degree relatives. Monitoring long-term outcomes in a larger 
number of first-degree relatives, combined with screening 
for genetic risk- factors such as bone morphogenetic protein 
type II receptor [BMPR2] mutations, could determine 
if these markers correlate with subsequent disease in 
susceptible populations. 

Several factors limit the interpretation of these results. While 
our study contained a relatively large cohort of Group I PAH 
patients, conclusions about the relevance of these markers 
to specific etiologies of PAH or to first-degree relatives of 
IPAH patients are limited by the relatively small sample 
size for each category. Our study exhibited a higher female 
predominance than other PAH cohorts such as the REVEAL 
registry (90% vs. 79%]. This relative enrichment of female 
patients may have been a function of the self-referral 
mechanism for enrollment, or due to chance. It should be 
noted that our PAH cohort is similar to other published 
populations of Group I PAH with regard to age, etiology 
of disease, functional status, and hemodynamic variables. 
Moreover, the control cohort was chosen to match the female 
predominance of the disease cohort (Table 1). In a small 
number of patients, we were able to measure most, but not 
all, of the biomarkers (NT-proBNP, sVEGFRl, and IL-15) due 
to insufficient quantities of blood. Therefore, correlations 
were obtained by comparing only patients for which both 
values were available (Table 4). Since the majority of 
patients in this study were enrolled in a nonclinical setting 
(Pulmonary Hypertension Association Scientific Sessions), 
clinical imaging, functional class, and hemodynamic 
assessments were obtained contemporaneously (±3 months 
by our definition) rather than simultaneously with 
blood sampling, and carried a significant missing data 
rate (Table 2). Given that CRP and OPG are elevated in 
other forms of vascular disease, further studies would be 
needed to control for the presence or absence of diabetes, 



cerebrovascular disease, and coronary artery disease as 
confounding variables. Our finding that sEng or sVEGFRl 
are useful diagnostic and prognostic tools for PAH will need 
to be confirmed in validation cohorts. 

In summary, this study has identified the potential role of 
two angiogenic modulatory proteins, sEng and sVEGFRl, as 
biomarkers of Group I PAH, and validated the role of CRP in 
the assessment and prognosis of PAH. Our findings on sEng 
and sVEGFRl provide further evidence of an important role 
for dysregulated angiogenesis in the pathophysiology of 
PAH. The ability of these novel biomarkers to detect PAH 
even in minimally symptomatic patients and to predict long- 
term survival suggests that they could reflect pulmonary 
vascular remodeling activity at earlier stages, and perhaps 
be used to expedite treatment during critical periods 
when the outcome of the disease may be most amenable 
to modification. 
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